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Abstract 

An explicit form of charged-lepton mass matrix, predicting m T = 1776.80 MeV from the 
experimental values of m e and (in good agreement with the experimental figure m T = 
1777.05^0 26 MeV), is applied to three neutrinos z/ e , z/,j, v t in order to correlate tentatively their 
masses and mixing parameters. It is suggested that for neutrinos the diagonal elements of the 
mass matrix are small versus its off-diagonal elements. Under such a conjecture, the neutrino 
masses, lepton Cabibbo — Kobayashi — Maskawa matrix and neutrino oscillation probabilities are 
calculated in the corresponding lowest perturbative order. Then, the nearly maximal mixing of 
Vn and v T is predicted in consistency with the observed deficit of atmospheric z/^'s. However, 
the predicted deficit of solar i/ e 's is much too small to explain the observed effect, what suggests 
the existence of (at least) one sort, vf \ of sterile neutrinos, whose mixing with v e would be 
responsible for the observed deficit. Perspectives for applying the same form of mass matrix 
to quarks are also outlined. Two independent predictions of |V^b|/|V^,| = 0.0753 ± 0.0032 and 
unitary angle 7 ~ 70° are deduced from the experimental values of | V us | and | V c b | (with the use 
of quark masses m s , m c and m&). In the last three Sections, the option of two sorts, ui and 

, of sterile neutrinos is considered. They may dominate neutrino mixing, and even cause 
that two extra neutrino mass states (arising then) are agents of some tiny neutrino instability 
and related damping of v e and Un oscillations. In Appendix, three conventional Majorana sterile 
neutrinos are discussed. 

PACS numbers: 12.15.Ff , 14.60.Pq , 12.15.Hh 
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1. Introduction 



In this paper, the explicit form of mass matrix invented for three generations of charged 
leptons e~ , \i~ , t~, and being surprisingly good for their masses [1], is applied to three 
generations of neutrinos u e) z/ M , z/ T , in order to correlate tentatively their masses and 
mixing parameters. This form reads 
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where the label / = v, e denotes neutrinos and charged leptons, respectively, while fi^\ 
£ (/) ) q,(/) anc i ^,(/) are rea i constants to be determined from the present and future 
experimental data for lepton masses and mixing parameters (/^ and are mass- 
dimensional). In our approach, neutrinos are assumed to carry pure Dirac masses. 

Here, the form (1) of mass matrices (M^if) an d {^aj) ma J De considered as a detailed 
ansatz to be compared with the lepton data. However, in the past, we have presented 
an argument [2,1] in favour of the form (1), based on: (i) Kahler-like generalized Dirac 
equations (interacting with the Standard Model gauge bosons) whose a priori infinite 
series is necessarily reduced (in the case of fermions) to three Dirac equations, due to an 
intrinsic Pauli principle, and (ii) an ansatz for the fermion mass matrix, suggested by the 
above three-generation characteristics (i). 

In the case of charged leptons, assuming that the off-diagonal elements of the mass 
matrix (m^J^ can be treated as a small perturbation of its diagonal terms {i.e., that 
o;( e )///( e ) is small enough), we calculate in the lowest perturbative order [1] 
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when the experimental values of m e and m M [3] are used as inputs. In Eqs. (2), the 

o o ( e ) 

first terms are given as m T = 6(351m M — 136m e )/125, ^ = 29(9m M — 4m e )/320 and 
"*= 320m e / (9m^ — 4m e ), respectively. We can see that the predicted value of m r agrees 
very well with its experimental figure m° p = 1777. 05^26 MeV [3], even in the zero 
perturbative order. To estimate (a^ / pL^^j , we can take this experimental figure as 
another input, obtaining 

{^f = 0.0241H, (3) 

which value is not inconsistent with zero. Hence, a^ 2 = 180^ MeV 2 due to Eq. (2). 

For the unitary matrix (U^j, diagonalizing the charged-lepton mass matrix [m^J) 
according to the relation = diag(m e , m M , m T ), we get in the lowest per- 

turbative order 
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2. Neutrino masses and mixing parameters 

In the case of neutrinos, because of their expected tiny mass scale jfM, we will ten- 
tatively conjecture that the diagonal elements of the mass matrix (il^) can be treated 
as a small perturbation of its off-diagonal terms (i.e., that fi^/a^ is small enough). In 
addition, we put = i.e., Nfff = 0. Then, we calculate in the lowest perturbative 
order the following neutrino masses: 
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are relatively small by our perturbative conjecture, while 
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As seen from Eqs. (5), the actual perturbative parameters are not £ and \i Du t rather 
£/7 and %/7, what is confirmed later in Eqs. (9). Note that < 0, the minus sign 
being irrelevant in the relativistic case, where only m 2 U2 is measured (cf. Dirac equation): 
|m„ 2 | may be considered as a phenomeno logical mass of u 2 . 
Using Eqs. (5), we can write the formula 
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which will enable us to determine the product from the observed deficit of atmo- 



spheric neutrinos u^, if ^ — ► oscillations are really responsible for this effect. 



We calculate also the unitary matrix \U^{) diagonalizing the neutrino mass matrix 
l^p} according to the rela 
perturbative order we obtain 



(m^J^ according to the relation C/Mt M^U^ = diag(m jyi , m V2 , m„ 3 ). In the lowest 
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with x = (125/2106)f = f /16.848. 

Denoting by v a = v e , , z/ T and ^ = vi , , z/3 the flavor and mass neutrino fields, 
respectively, we have the unitary transformation 

"a = E( Vt ) -^ = E^>- 
* — * \ /a* ' — ' 

i i 

where the lepton counterpart (Vi a ) of the Cabibbo — Kobayashi — Maskawa matrix is given 
as V = C/Mt[/( e ) ~ C/Wt or 

the approximate equality being valid for negligible / /i^ when ~ 5p a due to Eq. 
(4). Of course, in Eqs. (9) we wrote a = 1,2,3 for simplicity. From Eq. (10), we get the 
unitary transformation \u a ) = J2i \vi)V ia , where \u a ) = z4|0) an d \ u i) = are flavor 

and mass neutrino statesQ. 

In the limit of ^ v > — > (implying £ — > and x — ^ 0), we obtain from Eqs. (10), (11) 
and (9) the following unperturbed mixing formulae for v\ , Vi , U3 : 
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These display the maximal mixing between z/2 and in all three cases and a smaller 
mixing of v 2 — z^exp (i(p^^ with V\ in the cases of v e and z/ T , giving a minor 

admixture to v e and a dominating admixture to v T (in z/ M there is no admixture of v\). 

*In place of Vj, — J2 a Via^a one might use the notation v' a = J2p Vap v p, analoguous to d' a = Vapdp 
customary in the case of quarks where V a p — X) 7 (P^ v ^) a ^yjg- 
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3. Neutrino oscillations 

Once knowing the elements Vi a of the lepton Cabibbo — Kobayashi — Maskawa matrix, 
we can calculate the probabilities of neutrino oscillations u a — > vp (in the vacuum) making 
use of the general formula 

P(v a - up) = \(u fi \u a (t))\ 2 -E V JP V ^V*pV ia e 2ix ^ , (13) 
where \v a {f)) — ex P( — iHt)\v a ) and 

= 1.26693 Am^L/E , AmUfflJ.-^ , (14) 

if Am];, L and E are measured in eV 2 , km and GeV, respectively, with L = t and 
E — \p\ (c = 1 = h) denoting the experimental baseline and neutrino energy. 

It is not difficult to show that for the mass matrix (^M^p^J, as it is given in Eq. (1), 
the quartic products of V^ Q 's in Eq. (13) are always real (for any phase f^), if only 
V ia = U a v }* {i.e., Uj£ = 8 pa). This implies that P[v a -> vp) = P(isp -> u a ). In general, 
the last relation is valid in the case of CP invariance which, under the CPT theorem, 
provides the time-reversal invariance. Because of the real values of quartic products of 
ViaS, the formula (13) can be rewritten as 

P(u a ^^)=8p a -4Y / V J pV; a V*pV ia sm 2 x Jl (15) 

i<j 

without the necessity of introducing phases of these products. 

With the lowest-order perturbative expressions (9) for V ia = U%}*, the formula (15) 
leads to the following forms of appearance oscillation probabilities: 
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as well as of survival oscillation probabilities : 
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Thus, we get P (v e — > z/ e ) + P (z/ e — > z/ M ) + P (z/ e — > i/ T ) = 1 and two other obvious sum- 
mation rules for probabilities. Among these probabilities, P (v^ — > z/ M ) displays (in the 
lowest perturbative order) maximal mixing between v 2 and z/ 3 . 
In the lowest perturbative order, 



ar 3 i - ^21 = 2:32 = 14 (— £ + x) (l.26693|Mg } | L/E) 
due to Eqs. (8) and (14). Hence, 
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sin 2 £31 = sin 2 x 2 i + £32 sin 2^21 + £32 sin 2x 2 ± 
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in experiments where x 32 <C tc/2. When in such cases the relation (23) is inserted into 
the formulae (16), (17) and (20), its x 32 and x 2 2 terms can be neglected in the lowest 
perturbative order. 

Note that the mass formulae (5) imply m 2 ^ <C m 2 2 < m 2 3 , where m^/m^ U3 = £ 2 /49 3 + 
0(£ 3 ) and m 2 2 /m 2 3 = 1 - (2/7) (48£/49 + x) + 0(£ 3 ). Thus, the inequality x 31 > x 21 > x 32 
holds in all neutrino oscillation experiments (with some given L and E). 

We have calculated the neutrino masses, lepton Cabibbo — Kobayashi — Maskawa ma- 
trix and neutrino oscillation probabilities also in the next to lowest perturbative order. 
Then, in Eqs. (5) the mass m Ul gets no quadratic correction, while m V2 and Tfiy^ are 
corrected by the terms 
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respectively. Among the derived oscillation formulae, Eq. (20), for instance, is extended 
to the form 
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= 1 - (l - 0.00514 £ 2 ) sin 2 x 32 - 0.000816 £ 2 (sin 2 x 21 + sin 2 x 31 ) (25) 

displaying nearly maximal mixing between v 2 and v 3 . 

In the case of Super-Kamiokande atmospheric neutrino experiment [4], if — > v T 
oscillations are responsible for the observed deficit of atmospheric i^'s, we have x atm = 
a^32 *C X21 ~ what implies that sin 2 £21 = sin 2 x$\ = 1/2 due to averaging over many 
oscillation lengths. Then, Eq. (25) leads to the following effective two-flavor oscillation 
formula: 

P K = 1 - (l - 0.00350 i 2 ) sin 2 x 32 , (26) 

if we assume in Eq. (25) that 0.000816£ 2 = 0.000816£ 2 (2 sin 2 x 32 ) effectively. Identi- 
fying the estimation (26) with the two-flavor formula fitted in the Super-Kamiokande 
experiment, we obtain the limits 

1 - 0.00350 £ 2 = sin 2 20 atm ~ 0.82 to 1, 

Am 2 2 = Am 2 tm ~ (0.5 to 6) x 10" 3 eV 2 . (27) 

Hence, £ ~ 7.17 to and 

t-j- = 0.00334£ ~ 0.0239 to , 
a {u) fi (u) = 0.483Am 2 2 ~ (0.241 to 2.90) x 10" 4 eV 2 , (28) 

where Eqs. (6) and (8) are used. For instance, with sin 2 26> atm ~ 0.999 and Am 2 tm ~ 
5 x 10~ 3 eV 2 , we get £ ~ 0.535 and 



„(") 

0.00178, ~ 2.41 x 10" 4 eV 2 , (29) 

what gives the estimation 

a {u) ~ 0.368 eV , //("> ~ 6.55 x 10~ 4 eV . (30) 

Note that £ < 1 for sin 2 20 atm > 0.9965. As was already mentioned, our actual perturba- 
tive parameters are not £ and x, but rather £/7 and %/7 = 0.0594£/7. 

Having estimated and /i^-*, we can calculate neutrino masses from Eqs. (5) with 
(6) and (7). Making use of the values (30) (valid for sin 2 20 atm ~ 0.999 and Am 2 tm ~ 
5 x 10~ 3 eV 2 ), we obtain 

m vi ~ 2.76 x 10~ 4 eV , m V2 ~ -1.71 x 10 _1 eV , m Ua ~ 1.85 x 10 _1 eV. (31) 

Because of the smallness of these masses, the neutrinos z/ 2 , u 3 are not likely to be 
responsible for the entire hot dark matter. 

In the case of solar neutrino experiments, all three popular fits [5] of the observed 
deficit of solar z/ e 's to an effective two-flavor oscillation formula require Am 2 ol <C Am 2 tm 
what implies Am 2 ol <C Am 2 2 <C Am^ < Am^, if — > v T oscillations are responsible for 
the deficit of atmospheric v^s. Then, x so \ <C x 32 <C ^21 ~ ^31, giving sin 2 £32 = sin 2 X21 = 
sin 2 £31 = 1/2 due to averaging over many oscillation lengths. In such a case, Eq. (19) 
leads to 

193 

P (is e -> u e ) = 1 - j-^ = 1 - 0.0402 = 0.960 , (32) 

predicting only a 4% deficit of solar z/ e 's, much too small to explain solar neutrino obser- 
vations. 

An intriguing situation arises in the case of formula (16) for P (v^ — > v e ), if — > v T 
oscillations really cause the bulk of deficit of atmospheric z/ M 's. Then, for a new x new = 
x m *C X21 5 x 3 i (with some new L and E) we may have sin 2 £21 = sin 2 £31 = 1/2 due to 
averaging over many oscillation lengths and so, infer from Eq. (16) that 

1 2 48 2 

P K -> v e ) = ^ sin 2 rr 32 + ^^£ 2 ~ 0.0204 sin 2 x 32 + 2.29 x 10~ 4 , (33) 
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where £ 2 ~ 0.286 (what is valid for sin 2 2# atm ~ 0.999 and Am 2 tm ~5x 10~ 3 eV 2 ). Such 
a predicted oscillation amplitude sin 2 26> ncw ~ 0.02 would lie in the range of sin 2 26*lsnd 
estimated in the positive (though still requiring confirmation) LSND accelerator experi- 
ment on — > v e oscillations [6]. However, the lower limit Am 2 SND £ 0.1 eV 2 reported by 
this experiment is by one order of magnitude larger than the Super-Kamiokande upper 
limit Am| 2 < 0.01 eV 2 . On the other hand, the small predicted oscillation amplitude 
sin 2 2# ncw ~ 0.02 would not be in conflict with the negative result of the CHOOZ long- 
baseline reactor experiment on u e — > P M oscillations [7]. 

In conclusion, our explicit model of lepton texture displays a number of important 
features, (i) It correlates correctly (with high precision) the tauon mass with electron 
and muon masses, (ii) It predicts (without parameters) the maximal mixing between 
muon and tauon neutrinos in the limit fi^ — > 0, consistent with the observed deficit of 
atmospheric v^s. (Hi) It fails to explain the observed deficit of solar z/ e 's. (iv) It predicts 
new — > u e oscillations with the amplitude consistent with LSND experiment, but with 
a phase corresponding to the mass squared difference at least one order of magnitude 
smaller. 

In the framework of our model, the point (Hi) may suggest that in Nature there exists 
(at least) one sort, of sterile neutrinos (blind to the Standard Model interactions), 
responsible for the observed deficit of solar z/ e 's through v e — > oscillations dominating 
the survival probability P{y e — > v e ) ~ 1 — P{v e — > v^f 1 ) [8]. In an extreme version of 
this picture, it might even happen that in Nature there would be two sorts, and 
of sterile neutrinos, where vjfi would replace v T in explaining the observed deficit 
of atmospheric u^s by means of — > oscillations that should dominate the survival 
probability Piy^ — > v^) ~ 1 — Piy^ — > v^) [9]. In this case, the constant for active 
neutrinos might be even zero (however, very small would be still allowed). Such a 
model is discussed in Sections 5 and 6. 

For the author of the present paper the idea of existence of two sorts of sterile neu- 
trinos is fairly appealing, since two such spin-1/2 fermions, blind to all Standard Model 
interactions, do follow (besides three standard families of active leptons and quarks) [8] 
from the argument (i) mentioned in Introduction, based on the Kahler-like generalized 
Dirac equations. Note in addition that the v e — > and — > oscillations caused 
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by appropriate mixings should be a natural consequence of the spontaneous breaking of 
electroweak SU(2) x [/(l) symmetry. 

In Section 7, a possibility is considered that two extra neutrino mass states, whose 
existence is implied by two sterile neutrinos and cause in the Standard Model 
framework some tiny neutrino instability and related damping of u e and oscillations. 

4. Perspectives for unification with quarks 

In this Section, we try to apply to quarks the form of mass matrix which was worked 
out above for leptons. To this end, we conjecture for three generations of up quarks u , c , t 
and down quarks d , s , b the mass matrices {M^^j and {M^J^, respectively, essentially 
of the form (1), where the label f = u, d denotes now up and down quarks. The only 
modification introduced is a new real constant added to in the element 
which now reads 



M 33 ) = |^( 624 + ^ )+C7(/) ) • ^) 



24/jW 
25-29 

Since for quarks the mass scales ^ and fx^ are expected to be even more important 
than the scale fi^ for charged leptons, we assume that the off-diagonal elements of mass 
matrices (M^^ and (m^^ can be considered as a small perturbation of their diagonal 
terms. Then, in the lowest perturbative order, we obtain the following mass formulae 

,.{u,d) / (u,d)\ 2 



29 ' 

lj(u,d) 4 



,,(u,d) r)A 

"** = V1 ( 624 + £M + CM ) +BM (^3t) • < 35 > 



where 



AM = ^1 36 M = »M 10800 m 

29 320 - 5eM ' 29 31696 + 54C(«> C + 29e( u ^ ' V ; 

In Eqs. (35), the relative smallness of perturbating terms is more pronounced due to extra 
factors. In our discussion, we will take for experimental quark masses the arithmetic means 
of their lower and upper limits quoted in the Review of Particle Physics [3] i.e., 
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and 



m u = 3.3 MeV , m c = 1.3 GeV , m t = 174 GeV 



(37) 



m d = 6 MeV , m s = 120 MeV , m b = 4.3 GeV . (38) 

Eliminating from the unperturbed terms in Eqs. (35) the constants jjS u ' d ^ and e^ u,d \ 
we derive the correlating formulae being counterparts of Eqs. (2) for charged leptons: 
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320 V ' * ; 320 V J \ u (u 



2 



^ ^ ^ + » . ,39, 



The unperturbed parts of these relations are: 

o (u,d) 

6 . , // 24 o(u,d) 

m *' 6 = 125 ( 351mc > s ~ 136m M) + C 

{ N o («,d) 

21.9 1 _ Ar ^ 24 
1.98 ) GeV + ^9"25 C ' 

A* = (9m c s - 4m u d ) = <^ ^ r _ ^ MeV , 

320 ' ' \ 95.7 J 

o(M) 29m M f 0.0904 1 

£ = TM = { 1.82 ) • ^ 

In the spirit of our perturbative approach, the "coupling" constant a^'^ can be put zero 
in all perturbing terms in Eqs. (35) and (39), except for the numerator of 

the factor (a {u ' d) / ^ d) f that now becomes (a {u > d) / Ji '^f. Then, and are 

replaced by 



o (u,d) o (u,d) 

o(M) A* 36 o(M) A* 10800 
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o (u,a) (u,d) 

Note that the first Eq. (35) can be rewritten identically as m U:d = I 1 £ /29 according 
to the third Eq. (40). 

We shall be able to return to the discussion of quark masses after the estimation 
of constants and a {d) is made. Then, we shall determine the parameters and 

, „ o («) o (d) 

C 1 - > (as well as their unperturbed parts C and C ) playing here an essential role in 
providing large values for m t and mj. 

At present, we find the unitary matrices (U^^ that diagonalize the mass matrices 
[M^lf^ according to the relations U^ u,d ^ M^ u ' d ^U^ u,d ^ = diag(m Uirf , m CyS , m tb ). In the 
lowest perturbative order, the result has the form (4) with the necessary replacement of 
labels: 



(e) — > (u) or (d) , /i — > c or s , r — > t or b , (42) 

respectively. 

Then, the elements V a p of the Cabibbo — Kobayashi — Maskawa matrix V = C/( u )t 
can be calculated with the use of Eqs. (42) in the lowest perturbative order. Six resulting 
off-diagonal elements are: 



Ks ~ Vcd -29{m s e m c 6 ) ' 

29 \m b m t J 29 m b 

V ub ~ 16V3» (w) « W c »M-)-hpW) 
841 m c m b 

ft = (43) 
841 m s m b 

where the indicated approximate steps were made due to the inequality m t ^> m b and/or 
under the assumption that a^/m c ^> /m b [cf. the conjecture (46)]. All three diagonal 
elements are real and positive in a good approximation: 

v ud ~i- l -\v us \ 2 , ~ 1 - i|vy 2 - \\v cb \ 2 ,v tb ~i- l -\v cb \ 2 . (44) 

In fact, in the lowest perturbative order, 
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4 / / \ , , A 180° 

arg ^ ^77 sin 93 w - y9 w ^ - arg V cs , arg V tb ~ , (45) 

841 m c m s v 7 7r 

what gives arg\4d = 0.88° = — argV cs , if the values (46), (49) and (52) are used. 

Taking as an input the experimental value \V c b\ = 0.0395 ± 0.0017 [3], we estimate 
from the second Eq. (43) that 

29 

a (d) ~ mb \V cb \ = (355 ± 15) MeV , (46) 
8y3 

where m b = 4.3 GeV. In order to estimate also oS u \ we will tentatively conjecture the 
approximate proportion 

a (u) . a {d) _ Q(u)2 . Q(d)2 =4 ^ 

to hold, where =2/3 and = — 1/3 are quark electric charges. Note that in the 
case of leptons we had «M : a (e) = 0.37 : {Vim x 10 6 ) = 2.8 x 10~ 8 for the central 
value of [cf. Eqs. (3) and (30)], what is consistent with the analogical approximate 
proportion 

« M :a^~Q^ 2 :Q^ 2 = 0, (48) 
where = and = —1 are lepton electric charges. Under the conjecture (47): 

a («) _ ( 1420 ± 60) MeV . (49) 
In this case, from the second and third Eq. (43) we obtain the prediction 

2 

\V ub \/\V cb \ ~ — ~ 0.0753 ± 0.0032 , (50) 

29 m r 
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where m c = 1.3 GeV. This is consistent with the experimental figure | V^,& | / j V^j, | = 0.08 ± 
0.02 [3]. 

Now, with the experimental value \V US \ = 0.2196 ± 0.0023 [3] as another input, we 
can calculate from the first Eq. (43) the phase difference (p( u > — ip( d \ In fact, taking the 
absolute value of this equation, we get 
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COS 



8 m R 



1 + 16 



841 / m r \ 2 



4 



"MS | 



= -0.0301 



(51) 



with m c = 1.3 GeV and m s = 120 MeV, if the proportion (47) is taken into account. 
Here, the central values of and \V US \ were used. Hence, 



(pM - (pW = 91.7° = -88.3° + 180° 



(52) 



so, this phase difference turns out to be near 90°. Then, calculating the argument of the 
first Eq. (43), we infer that 



tanfargKs-V^)" 4 — T 
v ' m c 1 

what gives 



sin I <£>M — (pW 



4:(m s /m c ) cos (</?(") — ipW) 



-0.365 , (53) 



argK s = -20.1° + v? W . (54) 

The results (52) and (54) together with the formula (43) enable us to evaluate the 
rephasing-invariant CP-violating phases 



wg(V: s V; b V ub ) = 20.1° - 88.3° = -68.2° 



(55) 



and 



arg(V; d V*V td ) = -20.1° , (56) 

which turn out to be near to -70° and -20°, respectively (they are invariant under quark 
rephasing equal for up and down quarks of the same generation). Note that the sum 
of arguments (55) and (56) is always equal to (p^ — ip( d ) — 180°. Carrying out quark 
rephasing (equal for up and down quarks of the same generation), where 

arg V us -> , arg V cb -> , arg V cd -> 180° , arg V ts -> 180° (57) 
and axgV u d, argV cs , argV^ remain unchanged, we conclude from Eqs. (55) and (56) that 
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arg Kb -> -68.2° , argF td -> -20.1° . (58) 

The sum of arguments (58) after rephasing (57) is always equal to ip^ — ip^ — 180°. 

Thus, in this quark phasing, we predict the following Cabibbo — Kobayashi — Maskawa 
matrix: 

/ 0.976 0.220 0.00297e- i68 - 2 ° \ 

(V a p) = -0.220 0.975 0.0395 . (59) 

V 0.00805 e ^ 20 - 10 -0.0395 0.999 / 

Here, only \V US \ and \V c b\ [and quark masses m s , m c , m b consistent with the mass 

matrices (il^) and (M$)] are our inputs, while all other matrix elements V Q/ g, partly 

induced by unitarity, are evaluated from the relations derived in this Section from the 

Hermitian mass matrices (M^) and (M^J^ [and the conjectured proportion (47)]. The 

independent predictions are \V ub \ and wgV ub . In Eq. (59), the small phases arising from 

Eqs. (45), argKd = 0.9° and axgV cs = -0.9°, are neglected (here, arg (V ud V cs V tb ) = 0). 

The above prediction of V a p implies the following values of Wolfenstein parameters 

[3]: 

A = 0.2196 , A = 0.819 , p = 0.127 , r? = 0.319 (60) 
and of unitary-triangle angles: 



7 = arctan - = — arg V u b = 68.2° , (5 = arctan = — arg V td = 20.1° . (61) 

P 1 ~ P 

The predicted large value of 7 follows the present experimental tendency. 

If instead of the central value \V US \ = 0.2196 we take as the input the range \V US \ = 

0.2173 to 0.2219, we obtain from Eq. (51) y? (u) - ip (d) = 89.8° to 93.6° (with \V cb \ = 

0.0395 giving = 355 MeV), what implies through Eq. (53) that arg V us — (p^ = 

-20.3° to -19.8°. Then, after rephasing (57), argK& = -69.9° to - 66.6° and 

arg V t d = —20.3° to —19.8°. In this case, the Wolfenstein parameters are A = 0.2173 to 

0.2219, A = 0.837 to 0.802, p = 0.119 to 0.135 and r] = 0.325 to 0.312 (here, A^P 2 + V 2 = 

\V ub \/\V cb \ = 0.0753 is fixed). Thus, 7 = - arg Kb = 69.9° to 66.6° and f3 = - &rgV td = 

20.3° to 19.8°. 
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In contrast, if the central value \V cd \ = 0.0395 (giving = 355 MeV) is replaced by 
the input of the range V cd = 0.0378 to 0.0412 (corresponding to a (d) = 340 to 370 MeV), 
we calculate from Eq. (51) that - y? (d) = 97.3° to 84.9° (with \V US \ = 0.2196), 
what leads to argV^ — = —19.3° to — 20.9°. Hence, after rephasing (57), axgV u b = 
-63.4° to -74.6° and argV^ = -19.3° to -20.9°. In this case, the Wolfenstein 
parameters take the values A = 0.2196, A = 0.784 to 0.854, p = 0.149 to 0.0951 and rj = 
0.298 to 0.345. Thus, 7 = - arg V ub = 63.4° to 74.6° and (5 = - axgV td = 19.3° to 20.9°. 
Here, \V ub \ = 0.00273 to 0.00323 and \V td \ = 0.00738 to 0.00874. 

Eventually, we may turn back to quark masses. From the third Eq. (35) we can 
evaluate 

C M = ^-m tb - 624 - e M> - ^ ( ^) ' , (62) 

^CM) 24 ' /i( M ' d ) 24 \/i( u ' d ) J y ' 

what, in the framework of our perturbative approach, gives 



, „ o (u,d) 29 25 29 /o (u,d) o (u,d)\ ( a (u,d) 
C (u ' d) = C +-7^—mt,b —\5A -9B 



2 



o (u,d) OA t ' b o (u,d) \ o (u,d) 

/i 320 /i v ' \\t 

29 f (u,d) o(u,d)\ ( a {u4) \ 

A +B ) I TTwT ' (63) 



o(u,d) 1 J I (u,d) 



where 



o(«,d) 29 25 (u,d) f 4339 1 f 4340 1 tBA . 

C =7(M)24 m ^- 624 - £ ={ 733.2 } = { 733 j- ^ 

With the central values of and as estimated in Eqs. (46) and (49) we find from 
Eqs. (41) 

o (u,d) ( a( u ' d) \ 2 { 7 39 1 o (u,d) ( a M \ 2 f 2 66 1 

A (tm> J = { lf 6 } MeV • B (totJ = { ass } MeV • < 65 > 



where 



o(M) / , „\ 2 

^ / a (u ' d) \ f 65.6 



MeV . (66) 



We calculate from Eqs. (63) with the use of values (65) that 
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rM = J 4339 + 5.25 1 f 4344 1 f 4340 1 

" \ 733.2 - 49.5 J \ 683.7 J \ 684 J ' 1 ; 

Similarly, from the second Eq. (39), making use of the values (65), we obtain 

« M - { } ™ - { % } ™ - { Z } -v ■ m 

o(M) {u,d) 

We can easily check that, with the values (40) for At and e and the value (64) for 

o (u,d) 

C determined as above from quark masses, the unperturbed parts of mass formulae 
(35) reproduce correctly these masses. In fact, 



o (u,d) , . 

V o(u,d) 3.3 
m ud = £ = s n } MeV , 



mc ' s 29 9 ' ~ J { 120 j 

o (u,d) 

A* 24 / o(M) o(M)\ f 174 1 , N 

= ^ 25 ( 624+£ +C )={ 4.3 } GeV ' < 69 > 

The same is true for the unperturbed part of the first correlating formula (39). The - 

here omitted — corrections to Eqs. (69), arising from all perturbing terms in the mass 

formulae (35) (including the corrections from S/j,^ u ' d \ 5e^ u,d ^ and SC^ u ' d ^), are relatively 

small, viz. 

^m=| J 2 l * j MeV, <5m c , s = j _^ | MeV , <5m t , 6 = j ™ j MeV , (70) 
respectively. 

We would like to stress that, in contrast to the case of charged leptons, where m T has 
been predicted from m e and m M , in the case of up and down quarks two extra parameters 
and appear necessarily to provide large masses m t and (much larger than 
m T ). They cause that m t {m h ) cannot be predicted from m u and m c {m d and m s ), till 
the new parameters are quantitatively understood. 

Note that a conjecture about and might lead to a prediction for quark masses 
and so, introduce changes in the "experimental" quark masses (37) and (38) accepted here. 
The same is true for a conjecture about and ip( d \ 
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For instance, the conjecture that the phase difference (p^ — ip^ is maximal, 



v (u) _ v (d) = g o ? ^ 

leads through the first equality in Eq. (51) to the condition 

fm s \ 2 841 / m s \ 2 , l9 . , 

i + ieM --rf-^J K. 2 = o (72) 



\m r J 4 VaW 



predicting for s quark the mass 



m s = 118.7MeV = 119 MeV (73) 

(with = 355 MeV), being only slightly lower than the value 120 MeV used previously. 
Here, m c and mj are kept equal to 1.3 and 4.3 GeV, respectively (also masses of u , d and 

o(<0 o(d) o(d) 

t quarks are not changed, while H , e and C change slightly). Then, from the first 
equality in Eq. (53) 

tan farg V us - ip {d) ) = -4 — = -0.365 , arg V us = -20.1° + ^ {d) . (74) 
After rephasing (57), this gives arg 14b + arg Vtd = — V 9 ^ — 180° = —90°, where 

argK& = -69.9° , arg V td = -20.1° (75) 

ie., practically —70° and —20°. All \V a p\ remain unchanged (with our inputs of \V US \ = 
0.2196 and \V c b\ = 0.0395), except for \V t d\ which changes slightly, becoming 

\V td \ = 0.00814 . (76) 

Thus, in the Cabibbo — Kobayashi — Maskawa matrix predicted in Eq. (59), only \ V t d\ and 
the phases (75) show some changes. The Wolfenstein parameters are 

= 0.118 , ry = 0.322 (77) 

and A and A unchanged (here, the sum p 2 + rf = 0.118 is also unchanged). Hence, 
7 + p = 90° and a = 180° - 7 - f3 = 90°, where 
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7 = arctan - = - arg V ub = 69.9° , (3 = arctan — ^— = - arg V td = 20.1° . (78) 
P 1 - P 

So, in the case of conjecture (71), the new restrictive relation 



P 

holds, implying the prediction 



V 1 P 2,2 /^n\ 

, P +V = P (79) 



\v td \/\v ub \ = \ r 2 ' T = - = 2 - 74 ' ( 8 °) 



1 - p) 2 + rf _rj 
p 2 + r/ 2 p 

due to the definition of p and 77 from V„& and V t d- It is in agreement with our figures for 
I I and \ Vub\. Then, the new relationship 



1 m c 



4 m s a( u ) 



a^ u 'm s p 



(81) 



follows for quark masses m c , m s and Wolfenstein parameters p, 77, in consequence of Eqs. 
(43) and the conjectured proportion (47). Both its sides are really equal for our values of 
m c , m s and p, rj. 

Thus, summarizing, we cannot predict quark masses without an additional knowledge 
or conjecture about the constants /j,( u ' d \ e^ d \ C^ u,d \ a^ u,d) and (in particular, 

the conjecture (71) predicting m s may be natural). However, we always describe them 
correctly. If we describe them jointly with quark mixing parameters, we obtain two 
independent predictions of \V U } } \ and 7 = — argV^ the whole Cabibbo — Kobayashi — 
Maskawa matrix is calculated from the inputs of |T4s| and of |T4&| [and of quark masses 
m s , m c and m b consistent with the mass matrices (M^^ and {m^J)]- 

Concluding this Section, we can claim that our leptonic form of mass matrix works 
also in a promising way for up and down quarks. But, it turns out that, in the framework 
of the leptonic form of mass matrix, the heaviest quarks, t and b, require an additional 
mechanism in order to produce the bulk of their masses (here, it is represented by the 
large constants and C^). Such a mechanism, however, intervenes into the process 
of quark mixing only through quark masses (practically m t and m&) and so, it does not 
modify for quarks the leptonic form of mixing mechanism. 
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5. A model of texture with two sterile neutrinos 



Assume that there are two sorts, and z/M, of sterile neutrinos (blind to all Standard 
Model interactions and so, interacting only gravitationally). Conjecture that their mixings 
with two active neutrinos v e and respectively, dominate all neutrino mixings. Thus, 
five flavor neutrino fields, v a = u e , v T , i/M, exist in this texture and mix according 
to a neutrino mass matrix M'"'. This can be assumed consistently in the following 5x5 
form: 



MM = (M$ 



( 













\ 























mS 


























V o 













/ 



with AfM = MM 



|MM| andMM 



(82) 



|M^| exp (i</?M^ f or q, < ^ where the 



diagonal elements M-M, M.M and are given in terms of //M an d eM as in Eq. (1) 
(with f = v). Here, we put M$ = = M 5 H and even M^ = = M 2 H , the latter 
implying o;M = due to Eq. (1) (with / = i/). With such a specific ansatz as (82), 
all neutrino mixings are caused by the existence of sterile neutrinos responsible for the 
off-diagonal matrix elements My and M^'. 

It is important to notice that, according to the useful formula for electric charge, 
Q = J 3 L + Y/2 with Y/2 = I* + (B - L) /2, sterile neutrinos can carry no lepton number, 
L — 0. This may be a reason for M,M = = M^. On the other hand, the off-diagonal 
matrix elements M^ and m|', if nonzero, violate the lepton number conservation. 

The mass matrix of the form (82) leads to the following masses corresponding to five 
mass neutrino fields i/j = v\ , v 2 , ^3 , ^4 , ^5: 



rn 



m 



V3 



m 



"2, "a 



M 



± 



33 5 



M 2 M 



± 



M^' 



M 2 M' 



+ |mM| 2 , 



imMi 2 



'25 



(83) 
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Note that in Eq. (82) we used for simplicity a = l,2,3,4,5, which convention, if used 
properly, does not introduce any serious confusion with 2 = 1,2,3,4,5. 

The corresponding 5x5 unitary matrix U^ u \ diagonalizing the neutrino mass matrix 
(82) according to the relation = di&g(m Vl , m V2 , m V3 , m U4 , m U5 ) , takes 

the form 



If" = (U$) 










1 



VT+F2 




x 



-tip 







where 





1 





x 



VT+x^' 



1 



VT+x 1 







1 







VTW 1 



(84) 



X 



Y = 







+ \ 




( m{? 




2|M<?| 


V2|m£ } | 


m Vh - 


4' 








\Mg>\ 


2\Mg\ 


V2|M 2 ?| 



(85) 



Note that always < X < 1 and < Y < 1. 

The flavor neutrino fields u a are connected to the mass neutrino fields z/j through the 
five-dimensional unitary transformation 



"a = y E(Vi) ai v i (86) 

i 

with (V - ^ . = (y)* a = Via> where V = (V ia ) denotes the lepton 5x5 counterpart of 
Cabibbo — Kobayashi — Maskawa matrix: 



(a, (3 = 1,2, 3) 







, (87) 

5 afS (a,/3 = 4,5) / 

where (U a e p (cx, j3 = 1,2,3)) is the charged-lepton diagonalizing unitary matrix given 
perturbatively in Eq. (4). If there cx^/fj,^ (jointly with its numerical coefficients) is 
neglected, then ~ (5 Q/ g) and so, we can put in Eq. (86) 
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' / .«=^' , ) 1 „=( £/w,£/( "l,-Kl,=^ , • (») 

In our model, uj£ are given as in Eq. (84). 

6. Neutrino oscillations and their possible damping 

Having once found the extended Cabibbo — Kobayashi — Maskawa matrix V, we can 
calculate the probabilities P{y a — > vp) of neutrino oscillations v a — > vp (in the vacuum) 
i.e., the probabilities of (vacuum) oscillations of the flavor neutrino states \u a ) — > \vp), 
where \u a ) = f£ |0) and 

k) = EkW« ( 89 ) 

i 

with |z/j) = z/||0). If allowing that, in general, not all mass neutrino states \vi) are 
absolutely stable, then 

\ Ui (t)) = e'^-^lui) = l^e-^-^* , (90) 

where = ^Jp 2 + m 2 Vi ~ |p| + m 2 J ./2\p\ and 7, = (\m Vi \/E) 7^ are neutrino energies 
and decay widths (with 7^ and i? ~ |pj denoting the neutrino decay widths at rest and 
neutrino beam energy, respectively). Thus, generally, we obtain for neutrinos (in the 
vacuum) the following damped oscillation formulae: 

= ^a + E^a^^Ie^-^V^^-l] . (91) 

They are analogues of the formulae for K° — > and — > if oscillations. Note that 
Eqs. (91) imply the probability sum rules in the nonunitarity form 



£PK^) = ]T|VyV 2 ^, (92) 

p i 
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in spite of the unitarity of V. Of course, the rhs of Eq. (92) is equal to 1, if all (here 
involved) 7$ are zero. In this case, the damping in Eqs. (91) disappears and they become 
the conventional neutrino oscillation formulae. The same is true for the next Eqs. (93). 

If the quartic products in Eqs. (91) are real (as it turns out to be in our case), we can 
rewrite these equations in the form 

~ E VjWJWta sin 2 (5LZ^< t) , (93) 

3>i \ 4 J 

where the first term is equal to 



E y^y^v iol 

Writing (Ej-E^t = Am 2 ji L/2E and (jj + ji)t = (\m V] \~if ] + \m Ui \~if ] )L/E with Am^ = 
m 2 y . — m 2 ,., E = \p\ and L = t, and then expressing the neutrino masses m Vi and rest 
widths 7j-°' 1 in eV, the experimental baseline L in km and the neutrino beam energy in 
GeV, we can insert 



-(7j+7i)* 



(94) 



Ej - Ei AmjtL _ 



2 E 



(\ I (°) 1 I I (°)\T 

[Jj + li)t -> 5.07 ^ — = y j + yi (95) 

in Eq. (91) and (93) (here, c = 1 = k)f\ 

From Eqs. (93) with Vi a = , we derive in the case of our form (84) of U%> the 

following damped oscillation formulae for active neutrinos v e , , v T (in the vacuum): 

^The insertion L = vt with v = \p\/E ~ c (c = 1) is called by Lipkin [10] the "right handwaving" which 
converts the "gedanken oscillation experiment" in time into the real oscillation experiment in space. In the 
first experiment, a flavor neutrino is created by a weak-interaction source (of size < L) in a momentum 
eigenstate \v a ,p) being a superposition of a few energy eigenstates \vi,Ei) (with Ei = *Jp 2 + m 2 ,.) 
describing mass neutrinos evolving in time. Inversely, in the second experiment, the flavor neutrino is 
emitted in an energy eigenstate \u a ,E) given as a superposition of a few momentum eigenstates \vi,pi) 
(with \pi I = y 1 E 2 — m 2 . ) describing mass neutrinos propagating in space (the requirement of coherence 
within this superposition leads to the condition | \pi\ — \pj \ | <C 1 /source size) . In the first case E^ — Ej ~ 
Am 2 j/2\p\, while in the second \pi\ — \pj\ ~ A.m\JlE. Here, E ~ c\p\ (c = 1). A "wrong handwaving" 
would be the insertion L = ViU with Uj = p/ Ei. 
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P(v e > Vjx) = = P(z^ -> z/ e ) , 
P(z/ e -> i/ r ) = = P(z/ T -> i/ e ) , 
P(z/ M -> i/ T ) = = P(i/ T -> , 

= ( 1 + y2 j -(iT^) sin 2 (x 5 -a; 2 ) e -(— ), 
P(z/ r ^z/ r ) = e" 2s/3 (96) 

and those where, beside v e , z/ M , z/ T , the sterile neutrinos z/j e ) , i/^ participate explicitly: 



p(^>) = p(^-f 4 ) y + (^) 2 sin2(X4 _ xl)e -(, + ,), 

P(v e ^„V>) = 0, 

^K^i e) ) = o, 

P („ T _> = o, 

P(^ T _> i/M) = 0. (97) 
The probabilities (96) and (97) satisfy the sum rules (92) which now read : 



e~ 2yi + X 2 e~ 2yA 
P(„ e - „ e ) + P(„ e - „M) = , 

e -2y-2 i Y 2 e~ 2ys 

= — ^ — . (98) 

Note that damping in our neutrino oscillation formulae decreases with growing neu- 
trino energy E, because yi = 5.07|m^. \^f^L/E decreases. Thus, the larger z/ a -neutrino 
energy is explored in z/ a -neutrino experiments, the smaller damping influence is exerted 
on P{y a — > v a ), provided not all (involved) 7* are zero. Of course, the effect of damping, 
if any, is expected to be very small. 
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7. A mechanism of negligible damping 

Now, we turn to the discussion of a possible mechanism of neutrino instability i.e., 
instability of mass neutrino states. To this end observe that the neutrino weak current 

J {V) " = WEl^eL + VjlrfVvL + WIl^rL , (99) 

though it is diagonal in the active neutrinos v e , , u T , is no longer diagonal in the mass 
neutrinos v\ , v 2 , is 3 , z/ 4 , v 5 , if the sterile neutrinos vj, e \ is^ really exist. In fact, inserting 
in Eq. (99) the unitary transformation (86), we obtain generally, beside Wl^vu, some 
nondiagonal products vjZ^VjL (i ^ j), since only three of five products I^z/y'VaL are 
originally present in Eq. (99). 

For instance, in the case of our form (84) of the unitary transformation (86) with 
V* a = gives 



v e = . 1 U -Xv±e llp{v) \ , 
v* = , 1 (u 2 -Yu^ iv) ) , 



= - 



v. 



00 = (Yu 2 e'^ + v b ) . (100) 



Thus, in our case, the neutrino weak current (93) transits into the form 



(101) 

Since in the Standard Model lagrangian this neutrino weak current is coupled to the Z 
boson [with the coupling constant —g/{2 cos 6w) = — e/ (2 sin W cos 9w)], some neutrino 
decays of the type V; L — > VjV k v\ with = (1,4) or (4,1) and (2,5) or (5,2), and 



1 



1 + X 2 

1 

+ 



1 + Y 2 



Wl^il + X 2 T7^^u 4L -X(l 
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with similar (k,l), are Z-mediated, so that they can be real processes if only \m u .\ > 
\ m uj \ + \ m u k \ + |ra„J (here, Pi denotes an antiparticle of vj). 

In the case of our neutrino mass spectrum (83), we get the inequalities m Vl > |m„ 4 |, 
m V2 > \m VB \ and m V2 > m Ul , where in the last relation we make use of > . 

Further, |m„ 5 | > m vi , \m U5 \ > \m V4 \, m V3 > m V2 and > \m V5 \, if Y — X > /\M^\, 
Y> X,Y < (M$-M$)/\M$\ and Y < M$ /\M$\, respectively Thus, for Y-X > 
M$I\M$\ and Y < (M$ - M$)/\M$\ all these inequalities hold. In this case, 
therefore, 

\m Vi \ < m Ul < \m U5 \ < m V2 < m U3 , (102) 

showing that then \m V4 \ is the lowest neutrino mass. 

We can see that for any virtual decay v 1 — > z/ 4 u k v% we get 

m Vl - \m VA \ - \m Vk \ - \m n \ < m Ul - \m VA \ - 2\m VA \ = 2M[? - \J 2 + A\m[?\ 2 

= Mi? - 2\M$\X > or <0, (103) 

depending on X < or > /2\M$\. This implies that, a priori, the decay width of 
v\ neutrino may be 71 7^ or 71 = 0, respectively. Since \m VA \ < m Ul , no virtual decay 
^4 — > v\ v k ui can be a real process, what leads to 74 = for z/ 4 neutrino. 
Similarly, for any virtual decay u 2 — > ^5 v k Pi, we obtain 

m V2 - \m v& I - \m Vh \ - \m n \ < m V2 - \m v& \ -2\m U4 \ = M$ + M$ - \j M^ 2 + 4|M^| 2 

= M$-2\M[l ] \X > (104) 

if X < M$/2\M$\, where M$ = (4/9)(80/e^ — 1)M^ with e M < 1 (cf. Eq. (1) with 
f = v). If true, this gives a nonzero decay width 72 7^ for v 2 neutrino. On the other 
hand, for v 5 neutrino 75 = 0, since \m U5 \ < m V2 . 

Anticipating that 71 = (or is extremely small) and putting 73 = 74 = 75 = 0, 
we obtain from Eqs. (96) and (97) the following neutrino oscillation formulae (possibly 
damped if 72 7^ 0): 
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(IX \ 2 
sin 2 (x 4 - xi) = 1 - P(u e - z/( £ )) , 

P(z/ T -> !/ T ) = 1 . (105) 
Here, 

Xl -x 4 = 2.53 1 14 ' 11 , x 2 -x 5 = 2.53 1 25 ' 22 . (106) 

From the neutrino mass spectrum (83) and the definitions (85) of X and Y, we can 
derive the useful equations expressing M$ and \My±\ through X and Am 2 4 , as well as 
M$ and |M 2 ( 5 } | through Y and Am 2 5 : 

/l _ Y 2 \ 1//2 ( Y 2 \ 1 I 2 

Mtf = \L^ Am ^ , | M £>| = ^-A-Am^J (107) 

as well as 

/i _ y2 \ V 2 / V 2 \ 

M% = ^i-AjAm^J , |m£>| = ^4^Am 2 5 J . (108) 
Further, writing 

1 > (y^) 2 - sin 2 20® , 1 > (A^,) 2 = sin 2 26><"> , (109) 

we obtain 



l>X = tanfl (e) , 1 > Y = tan0 M , (110) 

where < 26^ < vr/2 and < 20^ < vr/2. We can see from Eqs. (108) that for a 
fixed finite l-M^I we get Am| 5 — > as Y — > 1, excluding in this limit the corresponding 
neutrino oscillations. On the other hand, if we insist in an argument to keep Am| 5 fixed 
and nonzero as Y — > 1, we formally have (M^ | — > oo, implying m„ 2 — > oo and \m Vh | — > oo. 
(In both cases — > as Y — > 1.) Analogical conclusions follow from Eqs. (107) for 
(Mffl and Am 2 4 (and M^) when X -> 1. 
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The first Eq. (105) enables us to ascribe the observed deficit of solar z/ e 's to v e — > 
oscillations. In fact, we can determine our parameters and |M^| putting 



IX \ 2 

" = sin 2 2# solar ~ 0.75 , Amj 4 = Am 2 olar ~ 6.5 x lO" 11 eV 2 , (111) 



.1 +X 2 

if the global vacuum fit to solar data [5] is chosen. Then, due to Eqs. (110) and (107) 

X = tanfl solar ~ 1/V3 = 0.577 , M$ ~ 5.70 x 10~ 6 eV , \M$\ ~ 4.94 x 10" 6 eV . 

(112) 

Here, we can see that M$ /2\M$\ = (1 - X 2 )/2X ~ 1/^3 ~ X. Thus, the condition 
leading to 71 = is satisfied on the edge [cf. Eq. (103)]. At the same time, this shows that 
the condition M$/2\M$\ > X, providing 72 ^ in the second Eq. (105), is fulfilled 
comfortably [cf. Eq. (104)]. 

Damping in the second Eq. (105) complicates our discussion, though it is natural to 
expect that this formula allows us to ascribe the observed deficit of atmospheric v^s to 

— > oscillations. In fact, anticipating that damping in this case is tiny [cf. Eq. 
(119)], we may write exp(— y 2 ) — 1 — y 2 and, therefore, 



/ 2Y \ 2 ( 2Y 



^ - sin 2 (x 5 - x 2 ) 



, (113) 



,1 + F 2 , 

where the coefficient at y 2 in the correction 0(y 2 ) is almost compensated to zero. Thus, 
we can put approximately 



/ 2Y \ 

[TTy^J ~ Sin2 29atm ~ °' 82 to 1 ' Am 25 - Am ltm ~ (0-5 to 6) x 10- 3 eV 2 , 

(114) 

where the recent data from Super-Kamiokande atmospheric neutrino experiment [4] is 
applied. Here, we will put, for instance, sin 2 26 , atm ~ 0.999 and Am atm ~5x 10~ 3 eV 2 as 
in Section 3. Then, 



Y ~ tanfl atm ~ 0.969 , M 22 ~ 0.126 x 10 _1 eV , \M^ \ ~ 1.99 x 10 _1 eV (115) 
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due to Eqs (110) and (108). 

Making use of the estimations (112) and (115), we can evaluate and ^ from Eq. 

(1) (With / = !/), 



e<"> = ^ ^-l.eixio- 2 , 

on A/f (") 

= fffgL „ 1.03 x 10- 2 eV , (116) 



and then, the neutrino masses m vi , m V4 , m U2 , m U5 and m U3 from Eqs. (83), 



m Vl ~ 8.55 x 10~ 6 eV , m„ 4 ~ -2.85 x 10~ 6 eV , 

m V2 ~ 2.05xlO~ 1 eV, m^ 5 1.93 x 10 _1 eV (117) 



and 



m 



M 24 



^ " ~29~25 ( 624 + eM ) ~ 2 ' 12 X 10-1 eV • (U8) 
These masses satisfy consistently the inequalities (102) and reproduce the experimental 
values (101) and (114): Am 2 4 ~ 6.5 x lO" 11 eV 2 and Ara| ~ 5 x 10~ 3 eV 2 . 

Now, we can evaluate the total decay width at rest, 7^, for a mass neutrino v, L decaying 
through the Z-mediated processes i/j — > Vj u k Pi, where Wj = Ej + E k + E\ > rrij +m k + mi 
with m n = \m Un \. In the case of WI2, and m 2 — m 5 dominating over m k and 
(A;, Z = 1, 4), we obtain the approximate formula 

72 ° } = 1 19^ (iT^) 2 (m2 " ms)4 (m2 + 2m5) ' (119) 

where the total decay width 7^ is the sum of four partial decay widths for v 2 — ■> 1/5 v k v\ 
with (/c , I) = (1 , 4) , (4 , 1) , (1 , 1) , (4 , 4) which are proportional to 

Y \ 2 / X \ 2 / Y \ 2 f X \ 2 ( Y \ 2 ( 1 \ 2 / Y \ 2 ( X 2 \ 2 



1 + yVU+xV ' Vi + rWi + xv ' Vi + rVVi + xv ' \i + y 2 J\i + x 2 J ' 

respectively, the sum of these weights being equal to F 2 /(l + F 2 ) 2 . In this calculation, we 
used the Standard Model coupling of the neutrino weak current (101) to the Z boson [with 
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the coupling constant — g/(2cos6 w ), where G F /y/2 = g 2 /(8M w ) = g 2 /(8M Z cos6 w )}, 
and considered the situation when (p 2 — ps) <C Mf at the rest frame of decaying z/ 2 : 
p 2 = (m 2 , 0). In Eq. (119), the factor 1/4 at the front is a consequence of using the 
neutral weak current (rather than charged weak current), while Y 2 /(l + Y 2 ) 2 stems from 
mixing of active and sterile neutrinos. 

If Y, m 2 and m 5 are estimated as in Eqs. (115) and (117), then the formula (119) 
gives (with the Fermi constant G F = 1.17 x 1CT 5 GeV~ 2 ) the extremely small value 

7 f ~ 1(T 59 eV (120) 

corresponding to the enormous lifetime r 2 = l/7 2 °' ) ~ 10 43 sec (as eV^ 1 = 6.58xl0~ 16 sec). 
This implies for the Super-Kamiokande atmospheric experiment that y 2 = 5.07m 2 7 2 °' ) L/i? 
~ 10" 55 with m 2 7^ 0) ~ lO" 60 , L ~ 1.3 x 10 4 and E ~ 1 expressed in eV 2 , km and GeV, 
respectively. Thus, practically, y 2 = and so exp(— y 2 ) = 1. If m 2 = m U2 and m 5 = \m V5 \ 
grow by one order of magnitude (what is the case when sin2# atm rises to 0.9999 and so, 
Y to 0.990), then 7 2 °' ) becomes not larger than ~ 10~ 54 eV and r 2 not smaller than ~ 10 38 
sec. 

Concluding the last Section, we can say that damping in neutrino oscillation formulae 
can be completely neglected, unless there are other sources of neutrino instability [11], 
more effective than the Z-mediated decays ui — > Vj u k P t considered in this paper. The last 
decays appear in the Standard Model framework if, additionally, there are sterile neutrinos 
mixing with the active ones and so, breaking the elektroweak symmetry SU(2) x U(l). 
Our discussion shows that the neutrino decay widths 7$ are zero for i — 1, 3, 4, 5 and are 
completely negligible for % — 2. However, our damped oscillation formulae (93) [and their 
more specific versions given in Eqs. (96) and (97)] can work for any sort of potential 
neutrino instability. 
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Appendix: Majorana sterile neutrinos 

The flavor neutrinos, three active v e , , z/ r and two sterile v^f 1 , i/^ , considered in 
Sections 5, 6 and 7, lead to five mass neutrinos v\ , u 2 , ^3 , ^4 , z/ 5 having pure Dirac masses 
(also in previous Sections neutrinos had always pure Dirac masses). Now, assume that 
there are solely three active flavor neutrinos, but they possess the "Majorana" 2x2 mass 
matrices 

each consisting of one Dirac and two Majorana masses, and m^' R \ respectively [12]. 
The mass matrices (A.l) imply the following mass term in the lagrangian: 

-C mass = \^{F IF) , (A.2) 

where 

i>) = u aL + [v aL f , i/W = z/ ai? + (z/ qR ) c (a = e , /i , r) (A.3) 

are the Majorana flavor neutrinos, three active and three sterile built up of chiral 
fields v a L, (j'ql) = (^cO'r and z/ Q j?, (z^ a i?) c = (v a ) L involved already in the Dirac flavor 
neutrinos u a = u aL + u aR and antineutrinos v c a = (u aL ) c + {u aR ) c . These conventional 
Majorana sterile neutrinos v$ contain, therefore, no extra neutrino degrees of freedom, 
in contrast to our previous Dirac sterile neutrinos v'f'^ = v^l^ + ^sR^ involving extra 
chiral fields vfjf* and vf^. Of course, in contrast to the Dirac, the Majorana neutrinos 
mix (maximally) the lepton number L. 

In the case of "Majorana" mass matrices (A.l), the overall neutrino mass matrix takes 
the 6x6 form 

M<„ = (^M«) = («„ ( ™P )) . (A.4) 

In this "pure-Majorana" mass matrix there is no mixing between flavor neutrinos from 
three lepton families a — e , n , r. 

Diagonalizing the "pure-Majorana" mass matrix (A.4), we obtain the neutrino masses 
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corresponding to six Majorana mass neutrinos 



z/ = cos0 Q i/W-sm0 Q i/W , 

: sin^z/W +cos^z/W , (A.6) 



where 



m Q D) 1 A m( l '-m|f'\ 1 

COS = ~ —;= 1 - 



with 0q, ~ 7r/4 + (m^) — /Am^ ~ 7r/4 (m£ may be negative). Here, the approxi- 

mate equalities are valid in the case of ~ m^. If in addition ~ ~ m^, 
then Eqs. (A. 5) give ~ and m^ 7 ~ 2vn!^\ In contrast, if ~ <C 
they imply m 7,// ~ fm^ (this case is known as the pseudo-Dirac case). Note that in 
the case of ~ the mass neutrinos and v 1 * are in an obvious analogy to the 
mesons K L = pK° — qK° and K s = qK° + pK°, where q/p ~ 1 — 2e ~ 1 is a counterpart 
of our tan# Q ~ 1 - (m^ - m^) /2m^ ~ 1. 

Any model with ~ "m^, leading to the nearly maximal mixing i/' 

(^I*^ T z/ i s ' ) ) /v^2, is orthogonal to the popular see-saw model with <C <C 



(R) 



which gives v J a ~ and i/^ ~ In fact, in this case we get from Eqs. (A. 5) and 
(A.7) 

< ^ * , m'J c mf) + ^j W ^ m<?> (A.8) 

and 

1 fm {D ^\ 2 

cos9 a ~ l - - « ~ l , sin0 Q ~ ~ . (A.9) 

2 \ m a / 



32 



In both cases, however, we may have very small m^. Notice that the present experimental 
limit on the (still not observed) neutrinoless double (5 decay (violating the lepton number 
L) allows for m^) of the order of 1 eV or smaller in both cases of ~ and 

m( L) < . 

With the use of the neutrino mass matrix (A. 4) we get the "pure-Major ana" oscillation 
formulae 



P - = K^le-^l^)! 2 = ^ Q sin 2 2^sin 2 - xQ (A.10) 

and 



P - ^) = K^le-^l^)! 2 = <W - P - . (A- 11 ) 

where x 1 ^ 11 = 1.27(m I c ^ II ) 2 L/ E with m 1 ^ 11 , L and E expressed in eV, km and GeV, 
respectively. Here, sin 2 29 a ~ 1 if mffi ~ . 

For a form of neutrino mass matrix more general than the "pure-Majorana" form 
(A. 4), more general mass spectrum and mixing appear. A fairly general mixing may be 
given by the following anticipated formulae for Majorana mass neutrinos: 

rt" = E uLfvL'" = E uit I c °1° t) 7 sin « a (A.i2) 

a a I Sin a ~\~ COS 6 a 

(with i — 1, 2, 3 and a — e , [i , r). Here = (U a u }^j is a 3 x 3 family unitary matrix 
diagonalizing a 3 x 3 neutrino family mass matrix = (M^p) through the relation 

(l/MtMMl/M) = Sij-rrii. 

If the Majorana mixing angle 6 a is taken as a universal 6 (what certainly would be 
the case for 6 a = 45° corresponding to = m^), then the mixing (A. 12) follows from 
the 6x6 neutrino mass matrix 

M H = « } ) with M# = M# ( ^ ^ ) , (A.13) 

all entries , A^ and A^ being dimensionless. In fact, such a form leads to the 6x6 
unitary matrix 

&<") = (0*)) with &M = UH ( 21 ) (A.14) 
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which diagonalizes according to the relation 



where 



(u^M^) t =5,H *„) , (A.15) 



m, 



l' n ^ mi \^ n with A ' = =f ^ ^ J +A( D ) 2 ~ =f A (D) 

(A.16) 

(i = 1,2,3) are neutrino masses. The approximate equality in Eq. (A.16) is valid for 
\(L) ^ \( R ) m Note that the mass matrix (A. 13) is the direct product of two matrices (3x3 
and 2x2) containing separately the family and "Majorana" degrees of freedom. Thus, 
also the spectrum (A.16) is multiplicative. 

In the case of neutrino mass matrix (A. 13), the "pure-Major ana" oscillation formulae 
(A. 11) are extended to the form 



P K a) - 4 ) ) = \{^ ) \e- lHt \^ ) )\ 2 = <W - sin 2 20£|^| 2 |^| 2 sin 2 (of-*?) 

i 

-^U^tu^U^uit {cos^sin 2 (af-xj) + sin^sin 2 (x 1 / -a?) 

+ cos 2 #sin 2 # [sin 2 (x 1 / -xfj + sin 2 (xj-x 1 /)] } (A. 17) 



j>i 



which holds when the quartic products of matrix elements U^} are real. In Eqs. (A. 17), 
x 1 ' 11 = 1.27 (ml' n ) 2 L/E. Here, sin 2 29 ~ 1 and cos 2 9 ~ 1/2 ~ sin 2 9 if A (L) ~ A™. 

The neutrino family mass matrix = (^M^J^ may be assumed in the form (1) 

(with / = i/). Then, in the case of small f = M$ /\M$\ and x = M$ /\M$\, the 
family unitary matrix = (U^^j is given in Eqs. (9). In order to derive from the 
neutrino oscillation formulae (A. 17) explicit results, we put A^ = A^ (= A^ M ^. In this 
case, the neutrino mass matrix (A. 13) has the form 

M ( "> = « } ) with M# = Mg ( ^ ^) ) , (A.18) 

and the neutrino mass spectrum gives = rrii where mj = are 

determined as in Eqs. (5) implying m 3 ~ |m 2 | >mi (m 2 = — |m 2 |). 
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With this mass spectrum, the further discussion depends on the ratio of A^ M ) and \( D \ 
We will consider two cases: (i) = \^ or (ii) \( M ^ <C \^ (the pseudo-Dirac case). 
We derive from Eqs. (A. 17) and (9) the following neutrino oscillation formulae: in the 

case {i) 



P { z/ a) -> z> } ) = 1 sin 2 1.27 * 

V e e > 49 V E ) 2- 

P = l_ s i n 2 (\.27 ^ 



49 2 

!m 2 A( D ) 2 L\ 



* - ^) - 4^ (l.27 4 ^-;^ (C>2L ) (A- 19) 

or, in the case (m) 



r,^fai faA Z 48 ^ 2/ _4m 2 A (M) A (D) L \ 387 

P - tf>) = 1 - ^ sin 2 (l.27^ A(g)L ) - sin 2 (l.27^ = , 

n / fa) raU 1 • 2 A ^4(m 2 -m 2 )A( D ) 2 L\ 1 . 2 / n Jm^^A^M 
P M°) -> j/W ) = — sin 2 1.27-^ — 3 sin 2 1.27 , 

V^ e /49V E / 2 • 49 V E ) 

(A.20) 

where the L's are three different experimental baselines. In these equations, the negligible 
constant terms come out from terms containing sin 2 of large phases averaged over many 
oscillation lengths determined by the leading terms with sin 2 of small phases. The phases 
in Eqs. (A. 19) and (A.20) were calculated in both cases from the relations 



{my 1 ) 2 - {ml' 11 ) 2 = m 2 (A*"' =f A^) 2 - m 2 =f X^f , 

{m^ 1 ) 2 - {ml' 11 ) 2 = m 2 (\W ± \Wf - m* (\W T X {D) f , (A.21) 

working for A^ = A^ (= A^ M ^. Note that the second and third Eq. (A.20) are not of 
the two-flavor form, in contrast to the second and third Eq. (A. 19). 

Comparing two first oscillation formulae (A. 19) with the results of solar and atmo- 
spheric neutrino experiments [cf. Eqs (111) and (114)], respectively, we get 
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48 

— <-> sin 2 2# sol ~ 0.75 , 4m 2 A (D)2 <-> Am s 2 ol ~ 6.5 x 10 -11 eV 2 (A.22) 
49 



and 



1 <-> sin 2 20 atm ~ 0.82 to 1 , 4m 2 .A (D)2 <-> Am 2 tm ~ (0.5 to 6) x 10~ 3 eV 2 . (A.23) 
Hence, we obtain 



'" ' (3.61 to 1.04) x 10~ 4 (A.24) 



m 2 



and, due to Eqs. (5), 



f = (49) 3 / 2 ^ ~ (12.4 to 3.57) x 10" 2 , (A.25) 

while m\-m\ = 14[(48/49)f + xll^l 2 ~ (1-80 to 0.52)|M 1 ( 2 ) | 2 with x = f/16.848. 
This estimation confirms that £ = M^/lM^I and x = -^22 V I -^12 ^ I are small. 

In contrast to solar and atmospheric results, the LSND result (cf. Ref. [6]), say, 
sin 2 2#lsnd ~ 0.02 and Am1 SND ~ 0.05 eV 2 cannot be explained in the case (i), since in 
the third Eq. (A. 19) 

4(m 2 - m 2 )A (D)2 < 4m 2 A (D)2 ~ (0.5 to 6) x 10~ 3 eV 2 < Am^ SND (A.26) 

for the estimation (A.25) (m 2 . ~ m 2 ~ 49|M^|). 

In the case (it), however, one may try to compare the third Eq. (A. 20) with the LSND 
result getting, say, 



-J- <-> sin 2 2^ LS nd ~ 0.02 , (m 2 - m 2 )A (D)2 <-> Am^ SND ~ 0.05 eV 2 . (A.27) 
49 

If in the case (ii) the relation 4m 2 ,A( M * ) A < - D * ) <-> Am 2 tm analogical to (A.23) holds approxi- 
mately [c/. the second Eq. (A. 20)], the comparison with (A.27) gives 

4m 2 AW Am 2 , , /A x 

atm - (0.1 to 1.2) x 10" 1 (A.28) 



(ml - m2)A( D ) Am£ SND 
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and 



A (M) 




Ami 



atm 



(9.2 to 2.6) x 10 



3 A ™1 
Am?, 



atm 



(0.92 to 3.2) x 10 



,-4 



\( D ) 



Am? 



'LSND 



LSND 



(A.29) 



since 




(A.30) 



through Eqs. (5) (in making the estimation (A.29) the value (A. 25) was used, which holds 
also in the case (ii) if Am\\^\^ <-> Am s 2 ol ). Thus, for the value, (A.29) of A^/A^ 
the third Eq. (A. 20) might be consistent with the LSND result. 

In conclusion of this Appendix, we can say that a simple neutrino mass matrix (A. 13), 
operating with three neutrinos u e , , v T only and being multiplicative in "Majorana" 
and family degrees of freedom, is consistent in a natural way with solar and atmospheric 
neutrino experiments, but not with the LSND result (that still requires confirmation). 
Such a consistency of "Majorana" option does not differ much from that based on the 
neutrino mass matrix (82) including two Dirac sterile neutrinos v^f 1 and These 
conclusions were drawn with the use of our family mass matrix (1) (with / = is), where 
the dominance of its off-diagonal elements was conjectured. The opposite conjecture 
of dominance of its diagonal elements does not change our conclusions essentially. The 
nearly bimaximal mixing that appears in the v^f 1 — > v^f 1 and — > oscillation 
formulae (A. 19) is a consequence of maximal mixings of with and r ith v[?>, 
reflecting the equality A^ = A^ and so, not holding in the see-saw model corresponding 



When discussing the Majorana flavor neutrinos and (a = e, /i, r), one pre- 
sumes that the superpositions (A. 3) defining formally these objects are really coherent in 
processes of electroweak interactions which operate on lefthanded chiral fields u aL = v^l, 
ignoring their righthanded counterparts v a R = v^r- 

The Dirac part of mass term (A. 2) and the kinetic term J2 a Vail • dv a can be expressed 
by v a as well as and v*g\ viz. 



to AW « A( D ) < \( R \ 
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-4?aL = E m^V- a v a = £ ro f ) („«„« + i/Wi/W) (A.31) 
and, up to the full divergence id ■ J^o^al^a, 

£kin = £w-^ = ^E (^7 • 9^ o) + ^7 " . (A.32) 

Thus, the deciding role in the coherence question is played by the Majorana part of the 
mass term (A. 2), 



= o E { m « L) { V <±L) C V aL + VaL (Vol,) 1 



m 



(R) 



(A.33) 



which can be presented also in terms of Dirac superpositions v a = u a L + v a R and v c a = 
{v a L) c + (v aR ) c , but only if = m^K Hence, if mffi ^ (or even if mffi ~ 
only approximately), the coherence of Majorana superpositions and v$ seems to be 
physically preferred over the coherence of Dirac superposition i/ a . 
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